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ABSTRACT 

It  is  shown  how  to  combine  incidence  matrices,  which  admit  Hermite- 
Birkhoff  quadrature  formulas  of  Gaussian  type  for  any  positive  measure,  in 
such  a  way  that  the  resulting  matrix  also  admits  Gaussian  type  quadratures  for 
any  positive  measure.  Moreover,  the  uniqueness  property  and  the  extremal 
property  of  the  formulas  corresponding  to  the  submatrices  are  transferred  to 
the  formula  admitted  by  the  composed  matrix.. 
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SIGNIFICANCE  AND  EXPLANATION 


Methods  for  approximating  the  integral  of  a  function,  given  the  values  of 
the  function  and/or  some  of  its  derivatives  at  several  points  in  the  interval 
of  integration,  are  investigated.  The  integral  is  approximated  by  a  weighted 
sum  of  the  given  data — a  quadrature  formula. 

It  is  shown  that  for  a  wide  class  of  different  data  configurations,  there 
exist  appropriate  points  of  evaluation  and  weights  such  that  the  resulting 
quadrature  formula  is  exact  for  all  polynomials  of  the  maximal  possible 
degree.  The  well-known  Gaussian  quadrature  formulas  represent  the  particular 
case  in  which  only  function  values  (and  no  derivatives)  are  employed. 
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A  NFW  CLASS  OP  HERMITE-BIRKHOFF  OUADRATURE  FORMULAS  OP  GAUSSIAN  TYPE 

Vira  Dyn* 


1 .  Introduction 

Recently  results  concerning  the  existence  of  quadrature  formulas  of  Gaussian  type 

related  to  Hermite-Birkhoff  interpolation  problems  have  been  obtained  by  several  authors 

f 5 1 ,  (1),  [2],  [31.  Given  the  incidence  matrix  E  »  {e. . }m+l  "-1  with  entries  consisting 

il  i-0,j«0 

of  zeroes  and  ones,  a  Hermite-Birkhoff  quadrature  formula  of  Gaussian  type  (HB-GOF)  is 

defined  as  a  formula  of  the  form 

b 

M.1)  /  pdo  »  l  a^p^'lx^,  p  6  n^_i  , 


such  that  the  number  of  parameters  in  (1.1)  equals  the  dimension  of  n  (the  space  of 

n— 1 

polynomials  of  degree  <  n  -  1): 


m+1  n-1 

n  -  I  I 

i«0  i=0 


eij  +  m 


In  (1.1)  do  is  a  non-negative  measure  supported  on  more  than  m  points  in  (a,h), 

<*0 . Vn>esB  where 

S"  "  {Y  "  (v0 . vm+1,l*  '  *0  <  V1  <  <  *m  <  Vn  “  b) 

n-1  n-1 

and  J  e, ,  >0,  1  <  l  <  m,  £  e,  >  0,  i-0  or  m+1. 

j-0  j-0  13 

This  notion  of  HB-QOF  extends  the  classical  notion  of  Gaussian  quadrature  formulas 


u  m 

/  pdo  -  I  aiP(xi ),  p  «  II 
a  i-1  2m_’ 
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where  a<x^<*,*<xm<b  and  a^  >  0,  i  =  1,...,m,  and  also  the  notion  of  the 
multiple  nodes  GOF  [6],  (7]. 

The  problem  of  characterizing  incidence  matrices  admitting  HB-GQF  is  posed  in  [41 .  and 

the  following  necessary  condition  on  such  matrices  is  proved: 

Result  A.  Let  E  -  "-1  admits  a  quadrature  formula  exact  for  n  ,  for 

-  ij  i»0,j»0  n-1 

some  non-negative  measure  do  supported  on  more  than  m  points  of  (a,b),  and  let  r  be 
the  minimal  number  of  ones  which  must  he  added  to  E  to  obtain  a  matrix  E  -  (e^  } 
without  odd  sequences  in  rows  1 , # • • ,  m.  Then  E  is  a  Polya  matrix  (the  number  of  ones  in 
any  first  t  columns  exceeds  t  -  1),  and 

m+1  n-1 

(1.3)  n  <  l  l  eii  +  r  . 

i-0  j=0 

In  (21  the  existence  of  two  classes  of  HB-GOF  satisfying  (1.3)  with  equality,  has  been 
claimed.  The  first  class  is  related  to  incidence  matrices  with  rows  1,...,m  consisting 
of  odd  Hermite  sequences  and  even  non-Hermite  sequences.  The  existence  of  this  class  of 
HB-GQF  is  proved  in  (11. 

The  present  paper  is  concerned  with  the  proof  of  the  existence  of  the  second  class  of 
HB-GQF,  corresponding  to  Incidence  matrices  which  can  be  decomposed  vertically  into  several 
submatrices,  each  admitting  an  HB-GOF  with  equality  in  (1.3).  It  is  also  proved  here  that 
such  a  formula  is  unique  and/or  has  an  extremal  property  in  case  each  of  the  HB-GQF 
corresponding  to  the  submatrices  is  unique  and/or  has  an  extremal  property  of  the  following 
type: 

Definition  1:  The  HB-GQF  (1.1)  related  to  E  -  {e,,}'!'+!,  has  an  extremal 

-  ij  i«0o»0 

property,  if  for  any  f  satisfying  f*n'  >  0  on  (a,b1: 


b  b 

(1.4)  min  ( ( -1 )M  /  (f  -  p)do)  -  (-1)ufj  fdO  -  [  "i j f ^ ’ (xi > 1 

pGP(E,f)  a  a  eij*1 


n-1 


where  w  -  I  e  and 

j-0 


-2- 


P(K,f)  -  {p|p  6  nn_1#p<:J  )(vi)  -  -  1,  Y  6  Sm}  , 

with  E  =  (e.  . )  the  matrix  of  Result  A  related  to  E. 
il 

The  existence  of  another  interesting  class  of  HB-GQP  has  been  established  in  f 3 ) » 

These  formulas  are  related  to  incidence  matrices  of  pyramidal  types 

(a)  Each  interior  row,  i,  contains  one  sequence  only  which  is  of  odd  lenqth 
startinq  at  column  k^. 

(b)  There  exist  1  <  I  <  J  <  m  such  that 

ki  +  li  *  ki_,  *  V  1  ■  2, ...» I,  kt  -  0,  i  -  I . J,  kt  ♦  »4  >  ki+1  >  ki.  i  -  J . .  -  1 

(c)  The  two  rows  i»0,  i-mtl  have  non-zero  entries  in  arbitrary  positions 

k  >  k,  and  k  >  k_  respectively, 
i  in 

The  proof  of  the  existence  in  [3]  is  based  on  the  extremal  property  of  Definition  1, 
characterizinq  several  other  classes  of  HB-GQF,  in  particular  the  classical  GQF  and  the 
multiple  nodes  GQF  (61,  (71.  This  class  of  matrices,  as  well  as  the  class  in  (1],  can  be 
used  in  composinq  new  HB-GQF  of  the  type  discussed  in  Section  2. 


2.  Main  Results 


<»iven  an  incidence  matrix  E  »  £e.  .}?+l  "  1  with  n  ones,  and  a  set  of  points 

13  3*0 


X  fi 

Sa,  we  denote 

by  p(£,X,f)  the  interpolating  polynomial 

from  0  , 

n-1 

to 

f  at 

the 

data 

<E,X): 

(2.1) 

p(E,X,f  >*3  *(xA)  »f^)(Xi),  eAj 

*  1  • 

Let 

E  be 

a  Polya  matrix  with  even  sequences 

only  in  rows 

1 , , 

and 

n-1 

{(n>  > 

let 

«  -  1 
j-0 

Bm+1 ,  j 

>  0. 

Then  It  -  p(E,X,f ) 1 (-1 )U  >  0  on 

[a,b]  whenever 

0 

on 

[a,b] . 

Proof : 

It  is 

sufficient  to  consider  functions  satisfying  f(n)  >  0 

on 

[*,b]. 

The 

case 

f(n)  > 

0  on 

[a,b] 

is  obtained  by  taking  the  limit  of 

such  functions. 

Since  E  has  only  even  sequences  in  rows  1,...,m,  all  the  Rolle's  zeroes  of  the 
derivatives  of  g  «  f  -  p(E,X,f)  do  not  coincide  with  those  prescribed  by  E.  This 
together  with  the  Polya  conditions 


k-1  m+1  n-1  m+1 

(2.2)  l  l  e  >  k,  1  <  k  <  n  -  1,  l  l  e  «  n 

j-0  i-0  13  j-0  i-0  J 


imply  that  the  numbers  of  Rolle's  zeroes  of  g*k*  is  at  least  [  £  (  £  e  -  l)]. 

n-1  m  j-0  i-0  13 

Therefore  g*n_1'  has  at  least  one  (  l  l  e^  -  (n  -  1)  -  1)  zero,  which  is  the 

j-0  i-0  J 

maximal  possible  number  since  g *  0.  Thus  the  zeroes  of  g  and  its  derivatives  are 

only  the  Rolle's  extension  of  (E,X),  namely  those  prescribed  by  (E,X)  and  the 

corresponding  Rolle's  zeroes.  In  particular  all  the  interior  zeroes  of  g  are  of  even 
multiplicities,  and  therefore  g  does  not  change  sign  in  (a,b). 

To  determine  the  sign  of  g,  observe  that  if  em+1  k  “  0  then  the  right  most  zero 

fl.  \ 

of  g'  1  in  (a,b),  4^,  is  equal  or  greater  than  the  right  most  change  of  sign  of 
q(k+1)  in  (a,b),  tk+1.  Therefore  g*k'  (x)g*k,f1  *(x)  >0  in  (b  -  4^,b) ,  while  if 
g*k*(b)  -  0  («m+,,k  ■  1)  then  g*k' (x  )g*k+1 * (x)  <0  in  (b  -  e^^.b).  This  together 
with  g*n*  >  0  implies  that  (-U^gfx)  >  0  for  x  C  (a,b). 

As  a  direct  consequence  of  Lemma  1,  we  obtain: 


-4- 


Corollary  1  ••  Under  the  conditions  of  Lemma  1,  the  Peano  kernel  K(t,x)  5  K(t,x|E,X) 
for  the  error  in  the  interpolation  prescribed  by  (E,X): 

b 

(2.3)  f(x)  -  p( E, X. f ) (x )  -  /  f(n)(t)K(t,x|E,X)dt 

a 

satisfies 

(2.4)  (-1)WK(t,x)  >  0,  (t,x)  e  (a,b)2  . 

With  these  preliminary  results  the  existence  of  HB-GQF  corresponding  to  vertically 
decomposable  matrices  can  be  obtained. 

Theorem  1i  Let  E  =  (e. . }m+1  n_l  be  an  incidence  matrix  which  can  be  vertically 

-  l]  1=0,3 “0 

decomposed  in  *  submatrices,  each  consisting  of  rows  and  n^  columns  respectively: 

l 

(2.5)  E  =  E(n})  9  E(n2)  9  •••  9  E(nz),  l  *  n  , 

i-1 

such  that  all  the  non-zero  entries  of  row  i,  1  <  i  <  m,  belong  to  one  of  the  submatrices 
only.  If  each  submatrix  E(n^)  admits  a  HB-GQF  exact  for  )!„  _1 ,  for  any  positive 
measure,  and  if  each  E(n^)  satisfies  Result  A  with  equality  in  (1.3),  then  E  admits  a 
HB-GQF  exact  for  IIn  ,  for  any  positive  measure,  with  nodes  X  »  (xo' * '  *  ,xm-*-1  ^  6  ^ ' 
where 

(2.6)  a"  «  tx  =  (X<’>  U  X(2)  U  •••  U  x(t))|x(;i)  e  s  ■*,  1  <  j  <  4}  . 

Proof:  It  is  sufficient  to  consider  the  case  4=2.  Let  E  »  Efn^  •  Efn^)  be 

m. 

related  to  E  »  E(n,)  9  E(n2)  by  Result  A,  let  Y  6  S  be  the  nodes  of  the  HB-GQF  for 

m7 

E(n, )  and  the  measure  d<7,  and  let  ZSS  be  the  nodes  of  the  HB-GQF  for  E(n2)  and 
the  measure  ui(t)dt  with 

b 

(2.7)  ui(t)  =/  dO(x)K(  t,x  |  El^  ),  Y)  . 

a 

That  w(t)  is  of  constant  sign  follows  from  Corollary  1. 


-5- 


Consider  now  the  set  of  fundamental  polynomials  for  interpolation  at  (E, X)  with 
X  m  Y*  '  *  U  2^ 1 

(2.8)  ■  4ivV  \v  -  1  ■ 


and  for  any  p  6  II 


p(E,X,f)(x)  -  I  f (;i)(xi)qij(x) 


P<x)  -  I  P1^ *(xi)qi . (x)  , 


b  . .  .  b 

/  p(x)do(x)  «  l  p  3  (xA  )  /  q^ (x)do(x) 
a  ~  a  3 

eij-1 


Now  {q  .  |e  .  »  1,  j  <  n  }  is  a  basis  for  II  ,  since  q.  .  6  II  ,  vanishes  on 

ij  l]  i  ni_1  13  n2-1 

(E<n2),Z),  and  therefore  is  identically  zero.  Using  the  HB-G^F  with  nodes  Y  6  S  1  and 

coefficients  (a  |e.  .  «  1,  j  <  n,}  admitted  by  E(n. ),  which  is  exact  for  II  . ,  we 
l]  13  1  1  n  -1 


(2.12)  /  q^UJdotx)  -  'i  akvq'^  (yfc)  *  0  if  e^  -  0  and  e^  ■  1,  0  <  j  <  n.,  . 


Consider  next  the  set  {q  |e  ,  •  1,  n.  <  j  <  n},  consisting  of  n,  polynomials  in  II  , 

lJ  *]  1  c  n- 1 

vanishing  on  (E(n.|) ,¥) •  These  polynomials  can  be  represented  as 

b  (n  > 

(2.13)  <3i j  (X )  -  /  qiJ  (t)K(t.x  |  E(n^  ).Y)dt  . 


b  b  (n  )  b  b  (n  ) 

(2.14)  /  q  (x)dC(x)  »  /  dt  q  (t )/  K(  t,x)  | E(ni  ),  Y)dc(x )  -  /  q.  (t)a(t)dt 

a  a  J  a  a  ^ 
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But 


(l>1  ) 


q.  .  G  H  for  e^.  =  1,  n^  <  j  <  n,  and  therefore  the  application  of  the  HB-GQF 

13  n  _ —  1  j 

*  m5 

admitted  by  Etnj)  for  the  measure  w(t)dt,  with  nodes  Z  6  S  and  coefficients 

"  ’*  0<  j  <  n2}'  yieldS 

b  (n  ) 

(2.15)  /  q.  .  u(t)dt  -  l  b  q.(  .  ’  (z.  )  -  0  if  e, ,  »  0  and  e  .  =  1 ,  n  <  j  <  n. 

a  13  ekv=1  k'''_n1  13  k  13  13  ’ 

n^<v<n 


Combining  (2.11)  with  (2.12)  and  (2.15)  we  conclude  that 


D  D 

(2.16)  /  pdo  -  [  [/  q  (x)dO(x)]p(3)(x.),  P  6  nn_,  • 


e  ,  . =1  a 
il 


Using  the  same  arguments  it  is  easy  to  show  that 


( 1 )  u  • ••  u  x*43 


Theorem  2;  E  -  Efn,)  9  •••  •  Efn^)  admits  a  HB-GQF  at  x  =  X1  1 '  U  •••  u  X  for 
the  positive  measure  da,  if  and  only  if  E(n^)  admits  a  HB-GQF  at  X*i3  corresponding 
to  the  positive  measure  u^(t)dt  with 


(2.17)  W  <t)  =■  1/  K(  t,x  |  E(n, )  9  •••  9  E(ni_.,  >,X(  1 3  U  •••  U  X(l-1))da(x) 


A  direct  conclusion  of  Theorem  2  is 

Corollary  2i  The  HB-GQF  admitted  by  E  ■  E(n^)  9  •••  9  E(n,j)  is  unique  for  any 
positive  measure,  if  and  only  if  the  HB-GQF  admitted  by  E(n^),  1  <  i  <  t,  is  unique  for 
any  positive  measure. 

Another  property  of  the  HB-GQF  admitted  by  E(n^),  i  =  1 ,...,£,  which  is  transferred 
to  the  HB-GQF  admitted  by  Etn,)  9  •••  9  Etn^),  is  the  extremal  property.  Note  that  (1.4) 
can  be  formulated  also  as 


b  b 

min  ((-1)M  /  (f  -  p)daj  =.  (-1)W  /  (f  -  p(E,X,f)Jdo  . 
p6P(E,f)  a  a 
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Theorem  3i  Let  E  satisfy  the  conditions  of  Theorem  1.  If  for  any  positive  measure 


f  >]do 


(-1  )U  I  [f  -  p(E,Y  U  5<2). 


b  b 


(2.21 ) 


.  (n  ) 

(-1 )“  /  /  [£  -  p(E,YU  CU,,f)]  (t)K(t,x|E(nl  ),Y)dtdo(x) 

a  a 


b  -./o)  <ni  * 

(-1)  /  w(t)[£  -  p(E,Y  U  5  ,f]  (t)dt 


where  u>(t)  is  defined  by  (2.7),  and  satisfy  (-1)  w(t)  >  0  on  [a,b].  Applying  tiie 


extremal  property  of  the  HB-GQF,  admitted  by  E(n2)  for  (-1)  u(t)dt  at  the  nodes  Z, 


obtain  for  g2  =  f^  1  -  p  1  (E,Y  U  £*2\f),  satisfying  g2  J  *  f'n)  *  the 


<n2)  ,  ,(n) 


inequality 


(-1)  /  f(-D 


(2.22) 


U,  b  U  (n  )  (n  )  _ 

'  M(t)]{f  -  p  (E,Y  ,f))dt  > 

a 

li  b  u  (n  ) 


>  (-1)  2  /  ((-1)  Vt))(f  1  -q)(t)dt 


<ni  1 


(ni>  (",> 

where  q_  «  p(E(n  ),Z,f  )  6  H  .  Observing  that  [f  -  p(E,Y  U  Z,f)J  =  f  -  q2, 

2  n2~  1 

we  finally  derive  from  (2.20),  (2.21)  and  (2.22)  that 


u  b  ..  Vj  b  „  (n  )  U 

(-1)  /  (f  -  p(E,5,f)]dO  >  (-1)  *  $  [f  -  p(E,Y  U  Z,f))  (-1)  wdt 

a  a 

=  (-1 )“  J  [f  -  p(E,Y  u  Z,f )]do  , 
a 


indicating  the  optimal  property  of  the  HB-GQF  admitted  by  E  at  Y  U  Z. 

The  following  two  matrices  are  examples  of  matrices  which  admit  HB-GQF  in  view  of 
Theorem  1  and  the  results  in  (1)  and  (3): 
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1  1  0  0  0  0 
10  1  10  0 
0  0  0  0  0  0 
0  0  0  0  0  0 


10  10  0  0 
0  0  0  0  0  0 
10  0  110 
0  0  0  0  0  0 


0  1  0  0  0  0  0  0 

0  1  1  1  0  0  0  0 

E2  =  00000000 

1  0  0  0  0  0  0  0 

0  0  1  0  0  0  0  0 


0  10  0  0 
0  0  0  0  0 
1110  0 
0  0  0  0  0 
0  0  0  0  0 


The  matrix  E2  admits  a  j^-GpF  with  the  extremal  property  of  Theorem  3,  in  view  of  the 
results  in  [3],  [6]  concerning  the  two  submatrices  of  E2> 


-10- 
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